Evidence for mass zeros of the fermionic determinant in 
four-dimensional quantum electrodynamics 

M. P. Fry* 

School of Mathematics, University of Dublin, Trinity College, Dublin 2, Ireland 

(Dated: April 20, 2010) 

Abstract 

The Euclidean fermionic determinant in four-dimensional quantum electrodynamics is considered 
as a function of the fermionic mass for a class of 0(2) x 0(3) symmetric background gauge fields. 
These fields result in a determinant free of all cutoffs. Consider the one-loop effective action, the 
logarithm of the determinant, and subtract off the renormalization dependent second-order term. 
Suppose the small-mass behavior of this remainder is fully determined by the chiral anomaly. 
Then either the remainder vanishes at least once as the fermionic mass is varied in the interval 
< m < oo or it reduces to its fourth-order value in which case the new remainder, obtained after 
subtracting the fourth-order term, vanishes at least once. Which possibility is chosen depends on 
the sign of simple integrals involving the field strength tensor and its dual. 
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I. INTRODUCTION 



Within the standard model fermionic determinants are required for the calculation of 
every physical process. These determinants produce an effective functional measure for the 
gauge fields when the fermionic fields are integrated. They are the means by which virtual 
fermion loops are incorporated into a calculation. Without them color screening, quark 
fragmentation into hadrons and unitarity would be lost. In quantum electrodynamics the 
biggest barrier to understanding its nonperturbative structure is its fermionic determinant. 
These determinants are therefore fundamental. 

They are also hard to calculate and physicists by and large lost interest in them during the 
1980s. With the advent of large machines lattice QCD physicists are now starting to include 
the determinant in their calculations. Analytic results for QCD and QED determinants are 
very scarce, especially in four dimensions. Such results as they become available will serve 
as benchmarks for determinant algorithms, including the various lattice discretizations of 
the Dirac operator in use, and hence a means of reliably estimating computational error, a 
major problem in lattice QCD at present. 

Most analytic nonperturbative results obtained so far deal with the dependence of the 
determinant on the coupling constant. Little attention has been given to their dependence 
on the fermion's mass. One notable exception is the work of Dunne et al. [1], which gives 
a semi-analytic calculation of the QCD determinant's mass dependence in an instanton 
background. 

In two-dimensional Euclidean QED the author has shown that mass can have a profound 
effect on its determinant. Namely, for a large class of centrally symmetric, finite-range 
background gauge fields the growth of the determinant in the limit mR << 1 followed by 

|e$| >> 1 is 



where det denotes the determinant, R is the field strength's range, and $ is the background 
field's flux In the massless case, the Schwinger model, the determinant is quadratic in 
the field strength. 

The second example of the nontrivial mass dependence of det in Euclidean QED2 is the 
presence of mass zeros. Let det be written as 




(1) 
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Indet = Ha + lndet3, (2) 

where II2 is the second-order vacuum polarization graph and Indets is a technical term, 
defined in Sec. II, for the remainder after the conditionally convergent second-order term 
has been isolated and made gauge invariant by some regularization procedure. Then there is 
at least one real value of m at which Indets = when < je'l'l < 27r , subject to some mild 
restrictions on the field strength [3!] . There may be other mass zeros. Now recall Schwinger's 
result 1^] that Indets = when m = 0. For fields with $ = then it is also true that |5| 



lim In dets = 0; (3) 

m=0 

otherwise not. So the result is this: when < |e$| < 27r the zero in m of Indets moves up 
from m = to some finite value m > 0. Beyond |e$| > 27r we can say nothing definite yet. 

The obvious question to ask is whether there are mass zero(s) in the remainder term of 
Indet in QED4 , denoted by Indets . The background gauge fields A^{x) considered in two 
dimensions have a slow l/\x\ falloff resulting in a nonvanishing chiral anomaly $/27r. Here 
we will consider a large class of 0(2) x 0(3) symmetric background gauge fields that also 
have a l/\x\ falloff with a nonvanishing chiral anomaly. If the small- mass behavior of the 
remainder is fully determined by the chiral anomaly, as in two dimensions, then there are 
circumstances in which mass zeros are present in the remainder. The idea of the proof is 
extremely simple: show that for m — )■ the remainder is negative and that as m — )■ 00 it 
becomes positive. The demonstration that the chiral anomaly determines the small-mass 
behavior of the remainder turns out to be nontrivial, and we are not able to settle this 
matter here. Evidence is presented that it does, but it is not conclusive. 

At this point it may be asked why these mass zeros for a special class of background 
gauge fields are of interest. First and foremost they are a truly nonperturbative result for 
the exact QED4 determinant. As such, they would serve as a benchmark result that lattice 
theorists could aim to reproduce. As discussed in Sec. II, once the second- and fourth-order 
contributions to det^en are isolated the remainder of detren is determined by the distribution 
of its complex zeros in the coupling constant plane. Little is known about how these zeros 
distribute themselves. The presence of mass zeros in the remainder terms in Indet^en must 
place a strong constraint on their distribution which future work could deal with. 
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In Sec. II detren is defined and some of its properties are reviewed. Section III introduces 
the background gauge fields used in the calculation. Section IV is an introduction to the zero- 
mass limit of the remainder and some of the subtleties involved. Section V establishes that 
all of the square-integrable zero modes of the Dirac operator p have positive chirality. In 
addition it is necessary to know the scattering states and low-energy phase shifts associated 
with the background gauge field, and this is done in Sec. VI. Section VII gives an analysis 
of the low energy behavior of the exact negative chirality propagator and seeks to justify a 
particular approach to proving that the chiral anomaly is sufficient to describe the small- 
mass limit of the remainder. Section VIII demonstrates that the remainder can become 
positive as m — )■ oo. Finally, Sec. IX summarizes our conclusions. 



II. QED4 DETERMINANT 

We begin by reviewing some established results for the QED4 determinant 

as. By 

fermionic determinant we mean the ratio of determinants of the interacting and free Eu- 
clidean Dirac operators, det{p — + m)/ det(/* + m), defined by the renormalized deter- 
minant on R^, namely 

det^en = exp(n2 + U3 + n4)det5(l - e^/), (4) 

where 



Indets = Tr 



ln(l - eSA) + E ^"^"^^^ 



n 

71=1 



(5) 



and S = (f + m)-^n2,3,4 are the second, third and fourth-order contributions to the one- 
loop effective action defined by some consistent regularization procedure together with a 
charge renormalization subtraction in 112 • The regularization should also result in Ha = 
by C-invariance, and it should give a gauge-invariant result for 114 • The remainder, dets, 
after these subtractions is gauge invariant and has a well-defined power series expansion 
without regularization. The remainder Indets in (12]) is given by ([5]) with the restriction 
n = 1,2. 

The operator 5/ is a bounded operator on the Hilbert space L^(M^, y/k"^ + rn^d^k) 
for y4^G L"(]R^), in which case it belongs to the trace ideal C„, for n > A 
[Cn = {ir|Tr(Ktir)t < 00}] This includes the case when A^{x) falls off as l/\x\ 
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as |x| — )■ oo. As a result det^ is an entire function of the coupling e, and it can be rep- 
resented in terms of the discrete complex eigenvalues l/e„ of the non-Hermitian compact 
operator S'/ jlo| : 



det5(l - eS4) = Yl 



k 



(6) 



By C-invariance and the reality of dets these eigenvalues appear in quartets ±en, ±& 



or as imaginary pairs. Because det^en has no zeros for real e when m ^ and 
detren(e = 0) = 1, it is positive for real e. Because S"/ G C„, n > 4, it is of order 4. 
This means that for suitable positive constants A{e), K{e) and any complex value of e, 



Idetr 



4 was that in 



exp {K{e)\e\ ~^^) for any e > 0. The first paper to show that detren is of order 



12|. 



In the coordinate space representation of the operator S'(P)/(X), the propagator is given 



by 



S{x) 



m Ki{m\x\) 
'iQ + m) , , . 



47r 



(7) 



Here S is an analytic function of m throughout the complex m-plane cut along the negative 
real axis. A theorem of Gohberg and Krem [13j states that if G Ci and is analytic 

in /i in some region then so is det(l — v4(/i)). In our case S/^ G C4+e , requiring the four 



subtractions from the logarithm in (151). These subtractions are easily incorporated into 

nn 

Gohberg and Krem's proof for 5/ G Ci , provided use is made of the inequality [lOj, 114| . 



|det„(l + A)| < e 



r^\\A\ 



(8) 



for A G Cn, \\A\\^ = Tr|y4M|2, and is a constant. Therefore det5(l — 5'/) is infinitely 
different iable in m on the interval (0, oo). 



The regularization procedure used here is Schwinger's heat kernel representation 15l |: 

1 _o , 



In detr 



(7F)t\ 



x) + 



247r 



2 ^^l> 



x] 



-tm' 



F^,{k)\^ / dzz{l-z)\n 



z{l — z)k'^ + 



+n4 + lndet5(l - 5/). 



(9) 



Here e has been absorbed into A^ , D"^ = {P — AY , a^^ = [7^, 7,/]/22, 7^^ = —7^ , F^^ denotes 
the Fourier transform of F^^, , and m is the fermionic mass. A second-order on-shell charge 



renormalization subtraction has been incorporated. All terms appearing on the right-hand 
side of iQ follow from the heat kernel expression on the left-hand side. The requirement that 
E f] L"'(M'^) and certain differentiability conditions on introduced later are sufficient 

n>4 

to ensure that makes mathematical sense. 
In the representation 



(10) 





^0 






f CTk] 




7o = -M 






|,7.= | 


[-Ok j 


, 75 = 7o7i7273 




one gets 



H+ 
H_ 



D^ + laF= I " I, (11) 



where 

H± = {P-Af-(T-{B±^). (12) 

Denote the remainder in ([9]) after removing the renormalization dependent second-order 
term by 

7^ = ^4 + lndet5(l-5/). (13) 

It will be shown in Sec.V that all the zero modes of the Dirac operator are confined to the 
positive chirahty sector for the class of 0(2) x 0(3) symmetric background fields used to 
calculate det^en. Differentiating ([9]) with respect to rn? allows one to isolate After some 
rearrangement of terms there follows 



m 



= lm'TiUH^ + m')-'-{H_ + m')-^ 



1287re 







(14) 



where the spin traces are now over 2x2 matrices. Equation f lMj) expanded as a power 
series defines the conditionally convergent fourth-order term in TZ. This requires iterating 
the second term in ffT4l) four times using the operator identity 

ft 



e 



t(p-+v) _^-tp^ = _ / dse-(*-^)(^'+^Ve-^^', (15) 
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with V = -AP - PA + + a ■ (E-B). The result is 



m 



f d^k E(fc)-B(-A:) + E(-fc)-B(fc) 



167r2 J (27r)4 z{l - z)k'^ + m? 



2 d ^_ 
+ m ^-^114 



— m^Tr 



-(AA^AyAyAVA + all perms, of A^, V ) 
+/lA^/\A^/\A^/\ . 



(16) 



The second term in f[T6|) is the remainder after adding the second-order contribution from 
the second term in ( 11^ to the last term. This remainder would be canceled by the first 
term in ( fT6|) were it expanded in a power series. It, as well as , were calculated from the 
regulated expansion of the second term in using (fTSll . 

The quantity Ilj is obtained by adding all the fourth-order terms in the expansion. It is 
the contribution of if_ to the photon-photon scattering graph. Its structure is 

nj = nr^'^^ + nr^^-^^ (17) 



scalar 
4 



where n|™^"'' is the contribution to IIJ neglecting the cr ■ (B — E) term in V . That is, 11; 
is the contribution to 114 in scalar QED4 multiplied by 2. The factor 2, and not 4, is due to 
the factor of i in the definition (j9]) of the spinor determinant. The remainder in f|T7j) is the 
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contribution to from the cr • (B — E) term in V. These two terms are 



n 



-'^^^ - ^ f d'p d\ r r^-- 

- -8^ y (2^0^ (2^ (2^ X " X ' 



1-21-22 



X 



X 



[zik"^ + Z2p'^ + z^q^ - {zik + Z2P + z^qf + m^] ^ 
1(1 - 2^i)2(l - 2z3fF,,{p - k)F,,{k)F^p{-q)F^p{q - p) 
+1(1 - 2z, - 2^2)'(1 - 2Z2 - 2z;fF^,{p - k)F^,{-q)F^^{k)F^f,{q - p) 
+1(1 - 2^2)'(1 - 2^1 - 2z2 - 2z3fF^,{p - k)F^,{q - p)F^p{k)F^p{-q) 
-(1 - 2zi)(l - 2^2) (1 - 2^3) (1 - 221 - 2^2 - 2z3)F^p{p - k)F^^{q-p)F^s{-q)Fsa{k) 
+ (1 - 2zi - 2z2 - 223) (1 - 2^2) (1 - 2z2 - 2^3) (1 - 2zi - 2^2) 
xF^^(p - k)F0^{-q)F^s{k)Fsa{q - p) 

+ (1 - 2zi)(l - 223)(1 - 2^1 - 222)(1 - 2^2 - 2z3)F^p{p - k)Fp^{-q)F^s{q - p)Fsa{k) 
+ three additional field strength terms. . 

(18) 



and 



n 



(T-(B-E) 



1 



d'^k d'^p d^q 



dziZi 



1-21 



dzo 



1-21—22 



dz-i 



X 
X 

+ 



TT-J (27r)4 (27r)4 (2vr)4 7o j, j, 

[zik'^ + Z2p'^ + z^q"^ - {zik + Z2P + z^q)'^ + m^] ^ 

(B - E)(g - p) ■ (B - E)(p - A;) - 1 [(1 - 2zi - 2z2 - 2zsf + (1 - 2^2)'] 

xF,,{p - k)F^,{q - p)] (B - E)(A;) ■ (B - E)(-g) 



[B - E)(A;) . (B - E)(p - A;) - i [(1 - 2^3)' + (1 - 2^1)^] F,,{p - k)F,,{k) 
x(B-E)(-g)-(B-E)(g-p) 
- (B - E)(-g) . (B - E)(p - A;) + 1 [(1 - 2z, - 2z2f + (1 - 2^2 - 2^3)'] 

xF^,(-g)F,.,(p - A;)] (B - E)(A;) ■ (B - E)(g - p) 
+ three additional field strength termsj. (19) 
In ( fTSll the additional terms are gauge invariant expressions such as 

{-q)Ft,u{q - p)qaFai3{p - k)Ffi^{k)q^, 
and in (fT9|) they are of a similar form, such as 

(B - E)(-g) ■ (B - E)(g -p)g„F,^(p - k)Fp,{k)q,. 



The full 114 graph has been calculated by Karplus and Neuman {l^, although the above 
results do not appear in their paper. Results f|T8|) and f lT9|) require extensive, but straight- 
forward, calculation. As a check on 114™''^^ we can specialize to the case of constant B and 
E. Using F^pFfs^F^sFsa = 2(B2 + B'^f - 4(E ■ B)^, Hf'^''^'^ reduces to 



j^scalar ^ 



'-4 



7207r4m4 



(^|(B2 + e2)2-(B-E)2), (20) 



where is a Euclidean volume cutoff. This is Weisskopf 's [17| and Schwinger's [15|] constant 
field result for scalar QED4's fourth-order effective Lagrangian continued to Euclidean space 
modulo a factor of —2. The factor 2 was discussed above. The minus sign arises from the 
difference in statistics: we are calculating a contribution to the scalar QED4 effective action 
from the spinor effective action. 

Continuing with our discussion of fll6p . A_ in the last trace is the exact negative chirality 
propagator (x|(if_ + m'^)^^\y) and A is the scalar propagator (x|(P^ + m^)~^|y). The 
regulating exponentials have been removed as the terms in the trace are fifth order and 
higher, and so the implicit loop integral is unambiguous. We leave the discussion of A_ to 
Sees. VI and VII. 

The insertions clutter up the trace in the sense that in a perturbative expansion of the 
trace they can be neglected. This is because they form part of gauge invariant expressions 
whose structure is already determined by the AP + PA terms in V. We do not see any 
justification for neglecting such terms in a nonperturbative treatment of the trace, but 
nevertheless this remark should be kept in mind. 

In order to discuss the m = limit of f lT6|) we must be more specific about the background 
gauge field. 



III. BACKGROUND GAUGE FIELDS 



QED determinants in constant field backgrounds have volume divergences and so are 
not defined on non-compact manifolds. Instead one considers the associated effective La- 
grangians, which do make sense. In the simplest case of the Euclidean QED2 determinant 
there is just a constant magnetic field, and the volume divergence arises from the degeneracy 
of the Landau levels. In constant-field QED4, by making two rotations (a Lorentz boost plus 
a rotation in Minkowski space) the operator {P — A)"^ can be transformed into the sum of 
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two two-dimensional harmonic oscillator Hamiltonians, leading to a degeneracy factor that 
grows as a four- volume [jj]. The lesson is that constant fields have too much degeneracy to 
define the determinant on a non-compact manifold. 

We have found that 0(2) x 0(3) symmetric background fields allow a satisfactory defini- 
tion of the QED4 determinant and that they are sufficiently tractable to permit substantial 
analytic analysis. Such fields were first explicitly considered in QED4 by Adler Q • In 



this paper these fields take the form 



18|42Q| 



(21) 



where Mn^ is chosen to be antiself-dual and is given by 



M, 



-l\ 



(22) 



/ 



This field has an 0(2) x 0(3) invariance, subgroups present in the reduction of 0(4) to 
0(3) X 0(3). It is further assumed that a(r^) is smooth, well-behaved at the origin, and 
satisfies 



a(r2) 



, r> R, 



(23) 



where z/ is a dimensionless constant. Without loss of generality assume z/ > 0. 

The orbital angular momentum operators of the first and second 0(3) subgroups of 0(4) 

are 



-(1) 



(2) 



d 


d 


d 


dxk 


dxo 


OXm 


d 


d 


d 


dxQ 


OXk 





which satisfy 



Spq'iUjkL^k^ , p, q, — 1, 2. 



(24) 



(25) 



The spin angular momentum operators in the representation ( fTOl) are 



(1) 




0(2) 




(26) 
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The total angular momentum operator relative to the second subgroup, 



7(2) _ r (2) , q{2) 



(27) 



commutes with ^: 



0, A; = 1,2,3, 



(28) 



while ^ is invariant only with respect to rotations about the third axis of the first subgroup: 



-'3 ' 



0. 



(29) 



We adopt the conventions of 22] for the four-dimensional rotation matrices -D'^^^^: 

Ki + ml,. 



-^3 ^-^mim2(-^) 



im2 y^)i 



mim2 



(30) 



The -D^^^2 are normalized so that 



m\m2 



{x)D. 



27r2 



2\2h 



(31) 



2/1 + 1 

where is the surface element in four dimensions. Some properties of these matrices appear 
in Appendix A. 



Following 22] we construct eigenstates of J*-^-* ■ J^^-*, 73^'' (eigenvalues J ± |, M) and J*^^^ • 
J^'^\ Jg^'' (eigenvalues j,m). In the positive chirality sector these are 








X] 



(32) 



and in the negative chirality sector they are 

/ 



^7. 



3,m 



XI 







— {j — m + l)2D' 
^ {j + m + l)^D- 



X 



[X 



(33) 
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( j + m)iD' 



p-2 



X 



(34) 



Due to ([28D and ([29]), eigenstates of ^ = f - / are of the form |22| 



(35) 
(36) 



^EjMmi^) = fir {x)+g{r (x), 

where the superscripts on "ipEjMm denote chirahty and E is the energy eigenvalue. In the 
following we will write 4'EjMm two-component spinors. 
From = ^e^'^^'^F^^ and ^ it follows that 

-16a^ — 16r^aa', 



and 



F F 



3 4 '2 



(37) 



(38) 



where the prime denotes differentiation with respect to r^. From (123|) and (137|) the chiral 
anomaly is 



167r2 



X F,jyF, 



(39) 



provided lim r^a = 0. Note, as expected, that is not square-integrable. But this does 

r->-0 

not matter as far as the remainder TZ in f[T^ is concerned. Recall that it is only required that 
An E f] L"(M^), which it does here. Furthermore, because we have chosen on-shell charge 

n>4 

renormalization the behavior of F^,^ for small k in the first term on the right-hand side 
of is regulated by the vanishing logarithm as — )■ 0. So everything in is finite. 

If one wishes to deal with a negative anomaly then instead of choosing M^i, in (12T|) 
antiself-dual, require it to be self-dual. For example, let 



-1 



(40) 
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IV. ZERO MASS LIMIT OF 7^ : PRELIMINARIES 



We will now discuss in a preliminary way the limit of ( !T6|) as m — )■ 0. Consider the first 



term. A working definition of the chiral anomaly for p on non-compact manifolds is 23 1. 

lim m^Tr [{H+ + " {H- + m^'] = I d'x *F^,F^,. (41) 

Because the manifold is a non-compact Euclidean one the right-hand side of f HTl) need not 
be the difference between numbers — n_ of positive and negative chirality zero modes. 
The remainder, if any, is related to the zero-energy phase shifts associated with H± 23 |. 
More will be said about this at the end of Sec. V. If the remaining terms in ( fT6l) vanish in 
the m = limit then (1391) and (|4T1) indicate that TZ in (|T3l) behaves as 

7?. ~ — Inm^ + less singular in . (42) 

m-i-O 4 

Thus, 7^ would become negative as m — )■ 0. 

A necessary condition for the vanishing of the remaining terms is that there be no 
zero modes in the negative chirality sector. It will be shown in Sec. V that this is true for 
our choice of gauge fields. Otherwise, A_ in ( !T6ll would develop a simple pole at m = 
and f H2]) would contain more terms varying as In for m — t- 0. But this is not a sufficient 
condition for the remaining terms in fITB]) to vanish at m = 0. 

One can see already from the second term in (fT6|) some of the subtleties involved. If B(x) 
and E(a;) fall off as , as our fields do, without any particular symmetry constraint then 
their Fourier transforms will be such that B(A;), E(A;) behave as as — )■ 0. In this case 
the integral will have an infrared divergence even when m 7^ 0. But this does not happen 
due to the 0(2) x 0(3) symmetry of the gauge fields. 

To see this define 

F^Ak) = [ d'xe-''^F,,{x), (43) 

J\x\>R 



with 



Then 



\xT>R V '^^ ^ ■ ^ ^ 



M,,.h{kR) + ^'^"^"^"^,^""^"^'^ (Jo(A:i?) + 2.h{kR)) 



(45) 
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and 

B>{k) ■ E>{-k) = -^^^^{J^{kR)UkR) + JlikR)). (46) 

Thus, B>(A;) ■ E>(— /c) behaves as R^/k"^ instead of l/k"^ as — )• 0. For large fc, F^^{k), 
calculated from an integral like (j^3|) but with |x| < R, behaves as sin(A;i? — 3Tc/A)/k^^'^ for 
any reasonable behavior of a(r^) near r = 0, such as a ~ Cr^ with /3 > — i or — i as required 
in Sec. VIII. Therefore, the integral in ( 1T6|) is absolutely convergent in the ultraviolet and 
its small-mass limit varies as (Inm^)^, allowing us to conclude that the second term in ( I16p 
vanishes in the limit m = 0. 

Now consider the third term in f|T6l) . m'^dll^ /dm'^ . Referring to (17), (18), (19) and fj45|) . 
simple power counting of momenta suggests that the integrals defining 114'^"'"'' and 114 
have a logarithmic mass singularity of the form (Inm^)" with n>l. If so, then the m = 
limit of m^dU^ /dm'^ would be nonvanishing, thereby falsifying (H2l) . It is encouraging that 
there is no immediate infrared divergence for m = that has to be canceled by the symmetry 
of F^j, , as in the second-order term of f|T6|) . The confluence of singularities in E and B is 
no longer present in fourth order. Nor are they present in higher orders due to the result 
cited in Sec. II that dets is well-defined for any yl„ G fl L"(M^), which includes our fields. 

n>4 

The fact that power counting does not ensure finiteness in the m = limit of 11 J indicates 
that the symmetry properties of F^y{k) will be required to give a finite limit to the individual 
terms in (18) and (19). Because of this reliance on symmetry, theorems on mass singularities 
of Feynman amplitudes known to the author are inapplicable here. The analysis required 
to give a definitive answer one way or another is beyond the scope of this paper. All we 
are able to do here is to present evidence for a finite limit of HI as m — 0. We note that 

n 

Adler's stereographic mapping to the surface of a 5-dimensional unit hypersphere |12l 
cannot help here due to the slow 1/r falloff of the vector potential. 

Consider, for example, the fourth term in (18). Let p ^ p + k, q ^ —q so that the chain 
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of field strengths is put into the form 



Fap{p)Fp^{-k -p- q)F^s{q)Fsaik) 



\F^p{p)Fcp{-k -p- q)F^M)FAk) 
\Fap{p)F^p{k)F^,{-k -p- q)F,,{q) 
B(p)-E(g) + B(g)-E(p) 



+ 



X B(-A;-p-g) ■E(A;) +B(A;) ■E(-A;-p-g) 



(47) 



The l/k"^ behavior of F^^{k) arises from the jQ{kR) term in 0451) . Fixing on the most singular 
terms, the first term on the right-hand side of 0471) contributes 



F^MhM ^ ^= g ^4 4 [{q ■ kf - {qoks - qik^ + gs^i - q^k^f] • (48) 



To isolate the leading singularity when A;, g — )■ we neglect the denominator in the fourth 
term in (18) when integrating over the angles defining . Using 



we see that the leading singularity for small k and q cancels. The first term on the right- 
hand side of ( H7|) also contains the term Fai3{p)Fai3{—k — p — q). The case k,p,q ^ with 
p « k,q reduces to the case just considered when the angles defining p^ are integrated 
over. The same conclusions follow for the small k and p behavior of the second term in ( H7|) 
as well as the case k,p,q ^ with q << p, k. 

Finally, consider the third term on the right-hand side of ( 1471) . Referring to the result 
( H6|) . the singularity in B(— A; — p — q) ■ E(A;) and B(A;) ■ E(— A; — p — q) is /k"^ at g = 
and not l/fc^. 

The gauge invariant expressions on the right-hand side of (147|) occur in all the terms in 
(18) and (19), and so the above cancellations occur there too. There are three additional 
field strength terms in the integrands of (18) and (19) and these must also be considered. 

Given that the second and fourth-order terms in (16) vanish at m = then presumedly 
so will all higher order terms m^dH^ /dim?, . . . generated by expanding A_ in (16). Since 
the zero modes reside in the positive chirality propagator A+ this expansion may have some 
justification. However, the scattering states extend down to zero energy, and these may 




(49) 
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result in nonperturbative mass singularities induced by A_. This will be examined in Sees. 
VI and VII. 



V. ZERO MODES 



In the representation (fTOj) p has the supersymmetric structure 



D 
-D^ 

and hence positive chirality zero modes are square-integrable solutions of 



(50) 



(51) 



where all subscripts on ip'^ have been dropped. From fl32p and fl35p 



(x) 



v 



(j-M+i)^F-(j + M + i)iG 



{j + M + l)-2F+U-M+l)-2G 



Dl. 1 (x) 



By (ETJ), G provided 



POO 

/ drr^^+\F^ + G^)<oo. 
Jo 



Inserting (!52|) in (15T|) results in 



G' 



2j + l 



{] + 1)2 - M^F -MG\= 0, 



(52) 



(53) 



(54) 



r'F' + (2j + 1)F + 



ar 



2j + l 



MF+ ./(j + D^-M^G =0. 



(55) 



Here j = 0, 1, . . . and — j — | < M < j + 1. Equations fl5^ and fl55|) appear in [2l|, |22| 
in a different notation, although the authors are considering an entirely different problem. 
There are three cases to consider. 
Case 1: M = -j - i. Then 



^+ = v/2jTT/^i,,„(x)G(r2) 







(56) 
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with 

and so 

G{T^) = G{rl)e'^^4''''^'\ (58) 

Since a = u/P for r > R, 'ip'^ G for j = 0, |, ■■■,jmax, where j^ax is the largest value of j 
for which z/ > 2j + 2 is satisfied. 
Case 2: M = j + | . Then 

1 



^+ = -V2j + lDi^{x)G{r') I J , (59) 



and 



^ = -G. (60) 

Based on case 1 it is clear that G^L^ for any u > 0. 

Case 3: \M\ < j + ^ ■ We claim that there are no zero modes in this case. To show 
this let z = r"^ , 

A = (2j + 1)F, (61) 



r = -2MF - 2^1 (j + 1)2 - M^G, (62) 



in dSH) and (l55l). Then these become 



Assume a has a power series expansion about z = and let 



+ (2j + 1)A = IzaT, (63) 



7T-1 

z^ = {2M + lza)A. (64) 






r = 5^c„^". (65) 
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Then A and F have the expansions 



A 



Co 



ao 



-z + 



ai 



+ 



Mai 



4(j + l) V2(2j + 3) 4(j + l)(2j + 3) 



(66) 



Man 



Mai 



+ 



M^a^ 



+ 



2(2j + 3) 4(j + l)(2j+3) 16(j + l 



(67) 



It will now be shown that the solution A , F that is finite at r = does not converge fast 
enough to make ip^ G for |M| < j + \ with i})^ given by ( l52l) . Let t = In^; and 

F = 7e-(^'+5)*, 
A = 5e-(-''+5)*, 

a = ae"*. (68) 



Then fl63l) and (EH) become 



(^7 



^ + (j + i)5 = i«7, 



(j + l)7 = (2M+i«)5. 



(69) 



(70) 



These are the same equations appearing in Eq. (5.24) of 21j. Following their analysis, 
multiply ( l69l) by 5 , ( |70l) by 7 and subtract: 



i^(7'-^') = (j + i)(f + 5') + 2M7<5. 



(71) 



Since 7 = r^-'^-^F, 5 = r^-'^-'^A and F, A are finite at r = 0, 7 and 5 vanish at r = 0. From 
fl53|) . if il)'^ G then F, G ^ ,"-2i-2-e^ e > and hence 7, 5 ~ for r — > 00. Integrating 

n!|) therefore gives 

dr 



r L 



(j + i)(7' + <5') + 2M75 



0. 



(72) 



Since |M| < j + |, f l72|) is impossible for real e. Hence the assumption that ■?/'"'" G for 
|M| < j + i is false. 

We now turn to the negative chirality sector. From ( l33l) . ( I34l) and ( l36l) . 



1 



'?/' (x) 



(j + m) 2 /}■ 



(j — m) 2 D 




-(j — m + 1)2/}' 



(j + m + 1)2D-'^^ (x) 




(73) 
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and ip G provided 

poo 

/ drr^^^^ [p + [r'^gf] < oo. (74) 
Jo 

From fl50l) negative chirality zero modes are solutions of 

D^- = 0. (75) 

Substitution of ( j73l) in ( ffSl) results in 



2^1 J -M+ \f-\lj +M + \ (2rV + 4(j + 1)^7) + -^/ J - M + \af-\lj + M + Ir^ag = 0, 

(76) 



+ M + i/'+'^J - M + 1 (2rV + 4(j + 1)^7)-'^J + M + ia/-yj - M + Ir^a^/ = 0. 

(77) 

There are again three cases. 

Case 1: M = -J - I . From ([73]), 

From ([77D 

2rV + 4(j + l)(7-r2a^/ = 0, (79) 



whose solution by inspection is 



9{r') = 9{r'o) [^^J e^^^^S (80) 
By dZH) ^- e L2 only if 

/•oo 

io 

and therefore g is too singular at r = to be in . 
Case 2: M = j + i . From ([73]), 



-(j - m + 1)2 D^.'^j i(x) 

(^+^+l)^<>+i(^) 

and from fl76|) . 

2rV + 4(j + l)^ + r^a^ = 0. (83) 
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As in case 1 g is too singular at r = to be in . 

Case 3: \M\ < j + |. We will demonstrate that ip'^L"^ . Let z = , 



y(j + M+i)(/ + r^g) + ^{j-M + r){r'g - /), 



A = W(j + M+i)(/-r^^7) + W(j-M + l)(/ + r^^?). 



Then ([76]), ([77]) become 



(j + 1)2 - M2 A = (M + lza)r, 



(84) 



(85) 



(86) 



+ I J + i + \/ (j + - r = (M + i^a)A. 



Making the expansions (l65ll gives for M 7^ 0, 

/ 

A = 60 



2M2+J + 1- J(j + i)2_M2 
1 H . _ ao-2 + 0(;z^ 



4M(j + 1) 



r = 6n 



v 



(j + 1)2 - M2 



1)2 -M2 



(87) 



v 



2(J + 1) 



/ 



(88) 



(89) 



For M = 0, 



A-^l + ^^^ + 0(.3)) 



(90) 



r = &o ( 77^^ + TTTT^^' + o(^^; 



^4(j + l)~ ' 2(2j + 3)~ ' 7 • ^^^^ 
Having established that there is a solution of and (|57|) that is finite at r = we now 
show that this solution is not square-integrable. There is also a solution that is too singular 
at the origin to satisfy fTMj) : we therefore ignore it here. Let 

A = z-^~^^5, 



= A (j + i)2-M2. 



(92) 
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Then §7^ reduce to 



z^-X5 = {M + lazh, (93) 
+ A7 = (M + laz)S. (94) 



Multiply ([93]) by 5 , ([94]) by 7 and subtract: 



where 



and 



j^(^^-f) = -(f + 5^), (95) 
ar r 



52 _ ^2 _ 4 [A(/2 - r^g') - 2Mr^fg] r^^+\ (96) 



5' + ^' = 2i2j + l)[f + ir'gr]r^^^\ (97) 



From (!7i)) . if 'i/' G then /, r^gf ~ r ^■'^ , e > 0, in which case \im{6'^ — 7^) = 0. 

r— >oo r=oo 

Because A , T are finite at r = 0, 7 ,5 = 0{r'^^~^^) as r — )■ 0. Hence, integration of ([95]) using 
([97]) gives 



poo 

{j + 1)2 _ M2 / rfrr^^+^ [/2 + {r^gf] = 0. (98) 
Jo 



But this is impossible for \M\ < j ' + 1. Therefore, the assumption that ip G L is false. 



Summarizing, it has been shown that all L zero modes of p have positive chirality and 
that these only occur when M = — j — i and for values of j satis fyin g u > 2j + 2. 



These results raise an interesting problem. The main result of 23| is 

1 

_ = n+ - r2_ + - 

/ TT 



n+ - + - ^/.(/) [5+(0) - 5r(0)] , (99) 



where n± are the number of positive and negative chirality zero modes. Si (0) are the 
zero-energy scattering phase shifts for H± in ([TT]) . /i(Z) is a weight factor, and / are the 
quantum numbers required to specify the phase shifts discussed in Sec. VI. We have just 
shown that n_ = 0. Suppose u = 3. Due to the condition for a zero mode derived above 
only j = 0, M = — I ,m = are allowed. So = 1, and it must follow that 

^ = i + ^EMO['5/'(o)-<5r(o)]. (100) 

Verification of this and ([99]) here would take us too far afield. 

21 



VI. SCATTERING STATES 



Having established that there are no negative chirahty zero modes it cannot be concluded 
that the m = limit of the last term in flTB]) is zero. Equation fHT]) when combined with f l39p 
and fl99p caution against this. They demonstrate that a particular zero-mass limit receives 
contributions from the scattering states of in (ITT]) . There seems to be no alternative 
to actually calculating the low-energy scattering states of H_ before deciding whether the 
m = limit of the last term in ( IT6|) is zero. 

Because p is anti-Hermitian we look for eigenstates of the form 

pilj = ikij. (101) 
Decomposing ip into its positive and negative chirality components and using fl50p gives 

D^D^/j- = (103) 

To get the scattering states ip~ it is easier to calculate and then use D'^ip^ = —ikip". In 
the representation (fTOl) the Zeeman term, |crF , is diagonal in the positive chirality sector, 
and so DD'^ = has the form 



DD^ 



Hi 

2 

H_i 



where the subscripts on H denote the eigenvalues of S^^ in (12^ . In (15^ let 



2j + l 



-2i- 



2^2 w . — ■ 2 

and decompose ip'^ into its upper and lower components: 



(104) 



(105) 



2j + i / Dii-iJ^) \ n^''^ 

27r2 



3 ' 
T 2 



2j + l 



27r2 



T 2 



where £ • a; = 1. Substituting Eqs. (11060 in turn in (I102p gives 



+ 



+ [AM ±2)a + r'^a^ ±r'^ 

dr 



(106) 



(107) 
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Equation f ll07p has to be supplemented by appropriate boundary conditions. For r > R 
according to f l23l) a = vjr^ . Let 



P±\=r--f±. (108) 



Then fll07p becomes for r > i? 



whose general solution is a superposition of Hankel functions 

/± = a±H'^^\kr) + l3±Hf\kr), (110) 

with 

A = [(2j + If + 4Mi/ + ^ (111) 

Choosing a± ,(3± so that 

PE,MM^-) / J cos (fcr - |(2j + 1) + 5;,^^, (A:) - |) , (112) 
gives for r > R 

PEair) = y|(e<^-f(^^--^^)+^i(^))if«(A:r) 

+e-<^-f (fcr)) , (113) 

where E = k"^ and a denotes j, M, ±|. The superscript on 5^ is a reminder that these are 
positive chirality phase shifts. The solutions flllSp are to be joined to the solutions of fll07p 
for r < R. This will determine the phase shifts. Equation flll3p fixes the normalization so 
that 

POD 

/ drpEa{r)pE'air) = 5{E - E'). (114) 
Then ip'^-^ in ( 11060 have the overall normalization 

{i^t:p^i^l,p)=5pp,5{E-E'), (115) 

where (3 represents j, M , m, ±|. 

The calculation of the low-energy phase shifts is outlined in Appendix B. Define the 
energy-dependent part of 5J by 

A+(fc) = ^-|(2j + l) + e(fc), modvr, (116) 
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and denote the expansion in powers of k of the logarithmic derivative of the interior radial 
wave function at r = i? by 

'rdrPEa 



PEo 



-f^-{kRYV^ + 0{kRf. (117) 



R 



The coefficients 7q, , are defined in Appendix B. Then for \M\ 7^ j + | 

tan A+ = -Xf^^^^^ + ^ t^^^)'' ' ^^^^^ 

with A > 1 for z/ > 0. There are several special cases to consider. 

Case 1: M = j + | and hence A = 2j + 1 + . From (11061) the only phase shift in this 

case is 5^. . 1 ik) and 



tan Aj;. , , , = -^^..^ ( ^ - 2r^,^.,.) (^y'^' (l + O [(A;/^)^ (fci?)^^-]) . 



r2(i + A) \\ + 

(119) 

Case 2: M = — j — ^ and hence A = |2j + 1 — z/|. From (11061) the only phase shift in this 



case is _i _i . 

h J 2' 2 



Case 2.1: 2j < u < 2j + 1 with < A < 1, 

2J + 1-U J \ 2 J ^ IK J ,K J }J 

(120) 

Case 2.2: 2 j + 1 < < 2 j + 2 with < A < 1, 

tan A+_^._ 1 1 = tan vrA + O [(fci?)^-^^] , (i2i) 
When A = I + e, |e| << 1, (I12ip becomes 

tan A- . . l + 2e -, 0ikRr 

''"'-^'-^ 7re + 2(r^._^._i__i -1) {fY '^ + 0(e2,e(fci?)i-2s(fci?)3-2e) 

Case 2.3: < u < 2j with A > 1, 
Case 2.4: z/ > 2j + 2 with A > 1, 

-^^t.4.. - xf?(A) f 2r,!:^'::^ ) )^"^ ^ ^^^^^ ' (^^)^^"^]) • 



(124) 
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Case 2.5: u = 2j + 1 with A = 0, 

tan A+_^_ , _ , = -vr (l - 2r^,_,._ i i ) (^^) \l + O [{kRf ln{kR)] ) . (125) 
Case 2.6: = 2 j + 2 with A = 1, 

^^"^ ^'-^-i-^ " In (f ) +^E- r,,_,_|,_ 1 + O [(A:i?)2 \n{kR)] ' ^^^^^ 

where 7^; is Euler's constant 0.577. . . 
Case 2.7: u = 2j with A = 1, 

tan A+_^._, = -vr (i - 2r,,_,_.,_.) (l + O [(fci?)^]) . (127) 

Although it is not required for the analysis here there is a compact relation between 
the interior wave functions and the phase shifts that ought to be mentioned, namely for 

2- / -J-y<^)P,Ml{s)p^M,^.{s) = Sin (5;^, (fc) - (A:)) . (128) 

J 

This is easily obtained by going back to fll07p and noting that 

2 



)' = -(rVpiP-i- (129) 



For r > R the right-hand side of (11291) vanishes. The constant p'ip_i — p' lPi in the region 
r > R can be calculated using (I112p . Then integrating (I129P from to i? gives (I128p . It 
holds for all energies. 

We now proceed to get the negative chirality scattering states, in particular / and g in 

by calculating D'^ip'^^ = —ikip'i . This results in two orthogonal states since 
^2 




v27r^A;r 



1 Ij + M+1 I ~{j~rn + l)iD;+^_, (x) \ f d 2j + ^ 



V2^krl^ 2J + 1 I (.+- + 1)^<L.(^) V^^'"" ^ 



;i3i) 
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EjMm,' 



Ij + M 



2j + l 



(j + m) 2 



M,m~ 




lj-M + 



ij-m + iy2Di'tl_, 



2j + l 



1 T^J- 



(j + m + 1)2L' 



M,m- 
1 




pEjM,-i{r) 



2j+r 

+ ar I p 

r / EjM,~l 



(132) 



These states are normalized so that 



(133) 



where /9 represents j, M, m, ±i. Because there are no zero modes in the negative 
chirahty sector we expect that the scattering states ( 11311) and fll32p form a complete set: 



oo i+5 j 

E E E 

i=0 -I m=-j 

2 



(134) 



VII. A_ AT LOW ENERGY 



The exact negative chirality propagator is 



^2 _j_ 



(135) 



with ip^^ given by fll3ip . fll32p and a = jMm, ±| . Now suppose A_(x,x') is divided into 
its low and high energy parts by replacing the integral in (1135^ by + J^, with AR << 1. 
Then our objective is to show that the low energy propagator has only minor deviations from 
the free propagator. This turns out to be the case except when u = 2j + 2 which results in 
a benign logarithmic mass singularity. The high energy propagator poses no obstacle to the 
m = limit in f|T6l) and is well-defined due to the assumed regularity of at the origin. 
In order to proceed we replace the differential equation (I107P with the integral equation 



P±(r) = ^±A/^^2j+i(A;r) + |v^ [ dr'V? 



(136) 
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where p± represents PEjM,±^ i ^± cire constants to be determined and 

da 



V± = {AM ± 2)a + r^a^ ± r 



dr 



(137) 



By differentiating f ll36p it can be verified that f ll07p resuhs. To fix A± require that p± 
join smoothly to the outgoing wave solution f lllSp at r = i? with 5J replaced by its energy 
dependent part defined in flll6p . Then 



P±iR) 



{Jx{kR) cos At {k)-Yx{kR) sinA+(A;)) 



(138) 



together with fll36p at r = i? determine A±. 

An upper bound on p +(r) for < r < i? will now be obtained. Starting with (2.60), 
(12.134) and (12.136a) in 2J] deduce that for z > and for fixed values of j 



(2j + l)!(l + z)2i+f 



(139) 



,(1) . ^, . ^ i2 fl + z 



m;uz)\<CH\i— 




(140) 



nz \ z 

where the constants Cj , Ch depend on j. From these results it follows that for z > z' > 

z \ 2i+i 



J2j+liz')Y2j+l{z) - J2j+l{z)Y2j+ 

with C of order one. Now iterate (11361) and let 

\p^±\r)\ = rip^J^\r). 
From fll36p and fll4ip the n*^ iterate satisfies 



< Cz 



Since < r < i? and we assume kR << 1, then f ll42p gives 



p^±{r)\ = r^|J^'{r) < 



(0), 



r f kr\ 



2j+l 



2\2) (2j + l)!" 



(141) 



(142) 



(143) 



(144) 
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Thus and UMl give 

4"' < 



(... ^ (¥)" fl^rY'^ \A^\ 



r2 



U4 



2^(2j + 1)! V 2 / 



TT 



ds\V±(s) 



(145) 



By (HH 



|p±WI < 



±1 



(146) 



2(2j + l)! 

valid for < r < i?, « 1. 

It remains to estimate the constants A±. Suppose \M\ 7^ j + |- From fll36p . fll38p and 
ffnSD with kR«l obtain 

£drV^ (1 - (^)''^' + 0{kR')^ l4(r)(|fcr)-2^-V±(r). (147) 
Because /c^ is an analytic perturbation of p± in f ll07p make the expansion 



P±[r) 



I kR\ ^ 

^^^(po±(r)+P2±(r)A;2 + 0(A;^)) 



(148) 



for < r < i?. Then to leading order in k 



A, 



kRy-''-' 



1 - 



r(2j + 2) (7±-A-i) 



r(A + i) (^^ + A-i) 



(149) 



Substitution of fll49p and fll48p into (11360 gives an integral equation for po± (^) whose solution 
we do not require here. The main conclusion is that for kR << 1, \M\ ^ j + ^ 

\-2j-l 



A 



± 



kR 



(150) 



where N±{j, M^u) is /c— independent. The centrifugal barrier term in (I107P for large j 
(certainly j » v) will cause p± to approach the noninteracting solution .,J^J2j+i{r-k) . 
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Hence, for large j, A — )■ 2j + 1, 7± — )■ 2j + |, A^±(j, M, z/) — )■ 1 and so A± — > 1. Equations 
(11461) and fllSOp give the explicit upper bound 



for < r < i?, kR << 1, |M| 7^ j + ^. The fc— dependence of this bound is consistent with 

dnHD. 

When M = j + case 1, Sec. VI, only p+ is relevant and there is no change in its overall 
k dependence. When M = — j — | only p_ is relevant, and the largest modification of f llSip 
occurs in case 2.2 when 2j + l<u<2j + 2 with 0< A < 1. Repeating the above analysis 
gives the same result as fllSip except that the factor {kR/2)^ is replaced with {kR/2)~^ and 
N_ is replaced with a new constant iV_. Thus, for case 2.2 fll48p is replaced with 

P-ir) = (Po- + P2-ir)k' + Oik')) . (152) 

The remaining cases when M = — j — | result in less singular fc— factors than {kR)~^. 

Now it is evident that the overall /c— dependence of p±{r) is not changed by differentiating 
it with respect to r. Therefore, the leading small /c— dependence of the radial wave functions 
in (I13ip and (I132p remains (fc-R)^, A > 1 for M 7^ ~j ~ \- Because of the factor k~^ 
multiplying them the negative chirality wave functions i'^-j^^j^i fall off as {kR)^~^ as — )■ 
for M 7^ — j — |, < r < i?. This statement was verified by deriving integral equations for 
the radial wave functions in (I13ip and (I132p starting from (11360 and proceeding as above in 
the derivation of the bounds on p±. 

The case M = — j — | has to be handled with care because when 2j + 1 < < 2j + 2 
we have noted that p_ behaves as {kR)~^ as A; — )■ and hence one might naively conclude 
that behaves as (/ci?)^^^^ with < A < 1 when k 0. This would induce 

a non-integrable singularity in the chiral propagator (I135p . This does not happen for the 
following reason. It may be explicitly checked that po- in (I152p is given by 

po-(r) = Cr^^+fe-^'-o^^^'^^^); (153) 

that is, it is a regular solution of (I107P when M = — j — | and k = 0. Here C and are 
arbitrary constants. This together with (11520 show that the relevant radial wave function 
in (11320 satisfies 



- ( - + ar - ^ ) p_(r) = 0(/.i^)-^ (154) 
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Hence, ibr^.,^ i continues to vanish as A; — )■ for M = — ?' — ^ and < r < i?. Another 

' ^ EjMm, — -^ ■'I — — 

potential non-integrable singularity from p+ in the second term in fll3ip is averted by the 
vanishing of the Clebsch-Gordon coefficient at M = — j — In view of the foregoing it is 
clear that the detailed analysis here is necessary. 

It is now required to examine the low energy behavior of the radial wave functions in 
f ll3ip and (11321) when r > R. We choose to deal with the most troublesome cases in Sec. VI 
first, namely those arising when M = —j — ^. We need only consider p_ in this case. For 
cases 2.1 and 2.3 when X = 2j + 1 — u get from (!23|) and flll3p 



^(|: + ar-^)p_(r) = - ^{J2,^,-Akr) cos Al^_, _^ 



1^2j+2-.(A:r)sinA+ . (155) 



For cases 2.2 and 2.4 when X = u — 2j — 1 then 

i + ar - ^) p_(r) = ^ {j.-2,^2ikr) cos A , 

- n-2,-2(A;r) sin A+_^._, . (156) 



Cases 2.5 and 2.7 are obtained by setting = 2j + l and 2j, respectively, in (11551) : case 2.6 is 

obtained from (fT56D by setting z/ = 2 j + 2. Using (fT20|l -(fT27D and Yp{z) r(p)(^/2)-77r 

for z obtain the following results with a denoting E,j, —j — |, m, — | 



case 


2.1: 




= 0{kRY+\^ 


< A < 1, 


case 


2.2: 




= 0{kRf-^, 


< A < 1, 


case 


2.3: 




= 0{kRf+\_ 


A > 1, 


case 


2.4: 




= 0{kRf-\, 


A > 1, 


case 


2.5: 




= 0{kR), 




case 


2.6: 




= 0(1), 




case 


2.T: C 


= o{kRy. 





(157) 



For case 1, M = ? + only ih^. . . i is relevant and 
1 f d 2j + -\ pr f 

^ (^^ - «^ J ^+ " ~V 2 (^2.+2+.(A;r) cos A^.^._,i_i 

-F2,+2+.(A:r)sinA+ , ,) . (158) 
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Thus . . 1 



0{kRy^^, A > 1. Finally, when \M\ ^ j + ^ and therefore A > 1, 
0{kR)^~^ based on ffTT3|) . ffTTSj) and the form of the radial wave functions 
appearing in (113 ip and fll32p . All of these cases are for r > R. 

Now return to the last term in f ll6l) and the interacting propagator fll35p . As noted 
earlier the centrifugal barrier term in f ll07p together with the regularity assumptions made 
on a(r) will cause the large j » v contributions to A_(x, x') to approach those of the 
noninteracting propagator. Therefore we need only consider a finite range of j in the search 
for a possible mass singularity in A_ that would result in a non- vanishing remainder at 
m = 0. 

For r » R the radial wave functions in fll3ip and fll32p are seen from flll2p to behave 

as 

I f d 27 + i\ II 



K \dr J. I 




Tck 



sin (^kr - i7r(2j + 1) + 6+^^^^, - Itt) 



+ 



cos{kr) 
kr 



(159) 



with k » 1/r, and similarly for the other group of wave functions. Therefore the leading 
large-distance behavior of these wave functions is the same as in the noninteracting case 
except for phase shifts. 

We have shown that the low energy wave functions V'^ -^^^j.! = 0(A;-R)''^~"^' , or less, for 
kR << 1 in the region between r = and r > R. Since A^ = (2j + 1)^ + AMu + u"^, then 
for z/ > and M > 0, A > 2j + 1 and 4'EjMm,±^ = 0{kRY^~^^, e > 0, which cannot lead 
to mass singularities in the last term in ( IT6l) that would result in a non-vanishing limit at 
m = 0. Hence we restrict our discussion to the case M < 0, in particular, the extreme case 



M 



-J - 



2' 



In (11570 the largest deviation from the noninteracting case is case 2.6 when u = 2j + 2. 
Focus on this mode in A_. This mode first opens up when u = 2, the threshold value of 
u for the formation of the first square-integrable zero mode in the positive chirality sector 
according to the discussion following (1581) . From (I135p . the second term in (11321) . (I156p and 
(I126p one obtains for the worst case r,r' > R 



A 



[x,x ) 



M{x,x') r^' 

Jo k^ + Tn'^ 



IT 



Mkr)Mkr') - -{Yo{kr)Mkr') 



kR\ 



+Yo{kr')Jo{kr)) / In ( — J + —Yo{kr)Yo{kr') ^ 



kR 



+ Ra, 
(160) 
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with AR << 1. The fc— independent matrix M is obtained from the second term in (I132p in 
the calculation of 'ipEai^)'^ELi^')^ -^a is the contribution to Ai from the region k > A. 
The most singular term in m in fll60p occurs in the first integral, written as 



^ ^j Uo{kr)Mkr') ^^kUkr)Ukr') 

A 



/j,2 _|_ ^2 J ^ f^2 _|_ 

POO 

Jo(mr<)ifo(mr>) - / (■), (161) 



where r<(r>) denotes the lesser (larger) of r,r'. The last integral in fll6ip can be put into 
the remainder in fll60p . Then for m — )■ the most singular behavior in m of A^{x,x') 
occurs in the mode j = ^z/ — 1 which has only a logarithmic mass singularity when r,r' > R 

—lu—i Mix x') 
A_ ^ (xjx') = —2 — ln(mr>) + less singular in m. (162) 

In summary, a mode-by-mode analysis of the exact propagator in ( 11351) uncovers only 
minor deviations from the free propagator in the low energy domain. If IIJ is finite at 
m = so that limm^ dU^ / dm^ = and if the role of the symmetry of Fn^ at large distances 

m=0 

in reaching this conclusion is well-understood then it should be possible to generalize this 
fourth-order result to m'^ dU^ / dm'^ , etc., obtained by expanding A_ in (fT6|) in a power 
series. We have shown in this section that its expansion is justified, considering that no 
nonperturbative singularities are induced in A_ by the scattering states that would cause 
the m = limit of the last term in f ITB]) to be nonvanishing. At this stage we do not see any 
other feasible way of demonstrating the vanishing of the last term in f lT6|) at m = since 
the available evidence relies on explicit gauge invariance and the long-distance symmetry of 

^ flU- 



VIII. LARGE-MASS LIMIT OF 7^ 



The leading term in the asymptotic expansion of TZ in (fT3|) for large m can be calculated 
from the effective Lagrangian density for QED4 in a constant field background 



This is possible provided F^^, is assumed to be smooth enough so that a meaningful derivative 
expansion of TZ can be carried out. Just how smooth will be made more precise below. 

The photon-photon scattering graph in TZ has been thoroughly studied by Karplus and 
Neuman 16|. Using their results or the comprehensive review of the Heisenberg-Euler 
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effective Lagrangian by Dunne [26|] one gets for large mass 



-0 ( ^ 1 d^xFapFapF^yd'^ F^^, — j d'^xF^j^adxFaydxF^pFp^ 



(163) 



/ d'x{F,,F^,f + >- f d'x{ *F,yF,yf. (165) 

Jr<R ^ 4 Jr<R 



The sixth-order graph obtained from the expansion of Indets is of order 1/m®. We now seek 
the conditions under which the leading term in 7^ as m — ?■ cxd, the right-hand side of f ll63p . 
becomes positive. 

From ( ]163p positivity requires 

J dSiF.yF^yf >^-j d'x{ *F^,F^,f. (164) 

From ([37]) and ([23]), *F^yF^y = for r > R and so by ([38]), F^^F^^ = Su^r^ for r > R. 
Then (I164p gives 

327rV4 ^ 7 

'r<R 4 Jr<R 

Referring again to f l37)) and f[38]) there follows the positivity condition 

/ drV6a'2[/a'2 + 4r2aa' + a^] > ^, (166) 

where a prime continues to denote differentiation with respect to r^. It is evident from (I166p 
that one class of fields satisfying the positivity condition is characterized by a steep rise in a 
in the region r < R where it obtains a maximum before descending as u/r'^ in order to join 
smoothly with a(r^) at r = R. Such an a is illustrated in Fig. [1] The class of admissible 
fields may be larger than this. 

In order for the remainder term in the asymptotic expansion in ( 1163^ to be finite it is 
necessary that F^^ be twice differentiable. From ( [2T1) the most singular term in F^j^ contains 



terms like x^M^aXaO'' and hence the most singular term in d'^F^y is of the form r^XyM^aXac" 
Thus, the finiteness of / F'^Fd'^F requires 







^d rf^^y^^^ 
dr J dr^ 



< oo, (167) 



and so a(r^) must be at least three times differentiable. For ease of analysis we assumed in 
Sec. V that a(r^) was regular at the origin, but this is not necessary. Condition (I167p only 
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R 

r 



FIG. 1: Sketch of a(r^) versus radial distance for a class of gauge fields satisfying conditions p66p 
and (fT671) 

requires a ~ Cr^ with /3 > — Of course requiring An G f] L"-{M^) rules out (3 < 0. Any 

r— >-0 



n>4 



branches in a(r^) away from r = of the form 

a{r^) ~ C(r2-r2)-, (168) 



must have a > 5/4 according to (11671) . 

Now it may happen that a given a(r^) does not satisfy condition (I166p . This could mean 
that either there are no mass zeros in the remainder defined by ( fT3|) or that there are an even 
number of such zeros. This cannot be decided here. In our search for definite information we 
go back to ( |T3i) and deal only with In dets, treating the photon-photon graph as a subtraction 
like the second-order graph. If fH2|) is true then the Inm^ singularity is from In dets alone. 
Then if the leading term in In dets 's asymptotic expansion in powers of 1/m is positive it 
certainly has at least one mass zero in the interval < m < oo. 



The leadirig 



;erm in the expansion of In dets in powers of 1/m is the sixth-order graph 
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given by [15, 

+0 J d'xF'F'F,,d'F,,, J d'xF^,^,d'F,,^^ . . . F^,^}j , (169) 

and hence the positivity condition is 

j d'x [13( *F^,F^,f - 8{F^,F^,f] F^pF^p > 0. (170) 
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From (l37Il and (l23|l this becomes 

/ d^x [13( - 8iF,,,F^,Y] F^pF^p > 

Jr<R 



10247rV5 



i?8 



(171) 



Further use of (!37|) and (!23l) results in the final positivity condition 



< (172) 



The most singular terms in the remainder of the asymptotic expansion in fll69p will arise 
from those containing d^F. Following the above discussion these will be finite provided 



drr 



da\^ d^c 



dr J dr^ 



< oo, 



(173) 



This requires a ~ Cr^ with f3 > —1/3, at least, and any branch points in a of the form 

r->0 

f ll68p must have a > 7/6. 

A necessary condition for positivity can be easily derived from Holder's inequality, namely 



I d'x\fg\ < (^j d'xin^y (^j d'x\g\^ 



1 1 

- + - = 1, P,? > 1- 
P Q 

Then with summation over indices understood, / = F^, g = ( *FF)'^, p = 3, q 



(174) 



d^x{ *FFfF^ < ( / d^xiF 



2 

( I d^xl *FFf^' 



(175) 



and so 



J d^x [13( *FFfF^ - 8(F2)3] 



< ( / d^x{F^^^ 



13 / rf^xl -8 / t/^x(F2^3 



(176) 



Thus, it is necessary that 



d X 



> ( ^ ) I d'x {F,,F, 



3 

flUj ; 



(177) 



for (11701) to be satisfied. It may be seen by inspection of (11721) that one class of fields 
satisfying it are those with a(0) ~ Nu/R"^, N > 2 and more or less monotonically decaying 
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R 
r 

FIG. 2: Sketch of a(r^) versus radial distance for a class of gauge fields satisfying conditions (|172|) 
and (fT73l) 

to p/R^ aX r = R as sketched in Fig. |2l Such fields will not satisfy the the positivity 
condition (11661) . 

To summarize, when fl42p . (11661) and (I167p are satisfied the remainder TZ in (1131) has at 
least one zero as m varies over the interval < m < oo. When fH21) . (11721) and (11731) are 
satisfied, Indets has such a zero. In this case the entire function in ([6]) somehow manages 
to reduce to unity at the mass zero(s). 



IX. CONCLUSION 

By choosing 0(2) x 0(3) symmetric background gauge fields we were able to make some 
provisional nonperturbative statements about the behavior of the Euclidean fermionic de- 
terminant detren of QED4 as a function of the fermionic mass. This determinant has the 
form 

In detren = na + n4 + Indets. (178) 

The second-order term contains a charge renormalization subtraction. The remaining terms 
are denoted by TZ in (fT3|) . It was assumed that for r > R the radial profile function a(r^) in 
fl2T]) takes the form u/r"^ for r > R, together with some mild regularity assumptions for a(r^) 
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for r < R. With these assumptions detren is free of all cutoffs, including the second-order 
term if on-shell charge renormalization is used. Then we showed that if the mass singularity 
of 7?. as m — )■ is fully determined by the chiral anomaly, then TZ has at least one zero as m 
varies in the interval < m < oo , provided conditions (166) and fll67p are satisfied. If not, 
then provided (172) and f ll73p are satisfied, Indets has at least one such zero at which the 
entire function in ([6]) becomes unity for any fixed coupling e. Then In detren is dominated 
by 114 for |e| >> 1, which is consistent with detren being an entire function of order four as 
discussed in Sec. II. If there is a mass zero such that TZ vanishes then In detren = II2 at this 
zero. If the number of mass zeros in TZ or In dets is even then they will not show up in the 
analysis here. 

This raises an interesting possibility. If << 1 then m does not have to be very large 
to make a meaningful l/m asymptotic expansion. So, presumably, there are one or more 
"small" mass zeros in the weak coupling domain. 

In plain language the result is this: set e^/hc = 1/137.... Select a gauge field that 
satisfies (166), fll67p or (172), fll73p . Adjust m until a mass zero appears. If m is the 
physical fermion mass then it probably does not coincide with a mass zero. But if, for the 
selected gauge field, m is near a zero then we would expect the remainder TZ or In dets to 
be anomalously small compared to the sum of the first few graphs in their expansion. By 
continuity there should be a class of gauge fields for which the physical coupling and mass 
coincide exactly with a mass zero. 

In establishing these results we also demonstrated a vanishing theorem when the field 
strength tensor is not (anti-)self-dual, namely, that all of the square-integrable zero modes 



of the Dirac operator are o 
of Brown, Carlitz and Lee 



one chirality. This is a generalization of the vanishing theorem 



27| . It would be useful to have a general vanishing theorem and 



to understand the physical principles underlying it. 

In Sec. VIII it was assumed that the expansion of TZ in powers of l/m is truly an 
asymptotic one so that the remainder after the series is truncated is of the order of the first 
neglected term. A proof is needed, but for the present it is an assumption physicists accept 
provided the background gauge field is smooth enough. 

Most of this paper deals with the question of whether it is indeed true that the leading 
mass singularity of TZ in flT3|) is determined by the chiral anomaly. We have presented 



evidence that it is. It is true for the case of constant B and E 



18| . but this is a formal 
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result as the determinant has to be made finite by a volume cutoff. And it is also true for 
the QCD4 determinant in the presence of an instanton background [28,]. It is evident that 
the analytic, nonperturbative analysis of four- dimensional fermionic determinants is still at 
an early stage and may yet yield some surprises. 
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APPENDIX A 



Here we list some properties of the four-dimensional rotation matrix. Several of the 
properties listed can be found in the Appendix of [22] from which our conventions and 
notations are taken. 

Let i = xo + ixs, T] = X2 + ixi. Then explicitly 



D^rrr.m.ix) = [{I - rm)Kl - m,)\il + m,)\il + m,)\]-^ J2 ^ ^ / , , , , (Al) 



— 'IIL2)'.{1 -t 'irtl)'.{l -t 'IH2)-\^ 7 

ni...n.4 

where = 0, 1, . . . and satisfy 

ris + Ui = I + mi, 
ns + rii = I + m2, 
Til + n2 + ns + = 2Z. (A2) 

The Dl^^j^^ are normalized according to ( 13T1) where we set = x^ + xf + x^ + x^. They 
satisfy 

A useful relation is 

I 

J2 D':.^^^{x)Dl^^^{x) = I. (A4) 

m2=—l 

The -Dmim2 satisfy an addition theorem 

Dl^^^{x)Dt,^^{y) = CUx-y), (A5) 

nii,m2=— / 

where the Gegenbauer polynomial is 

^1 / N sin(2/ + l)ip . ^ ^, 

sm(/? 

They also satisfy a completeness relation 

^ oo I 

— 5^(2/ + l) Dl^^2{x)Dt,^^{y) = 5{n,-ny). (A7) 

Z=0 m\m2=—l 
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The raising and lowering operators of the 0(3) subgroups give 



+ 1) - m2(m2 + 
V^(^ + l)-m2(m2-l)i^Li,m.-i- 

Other relations we found useful are 

^Jj-M+ \iD^^_.Jx) -^j+M+ \vD^^^. Jx) = v7T^r^D-^_,(a;), 
as well as 



r(l) r)' 
r(2) r)i 
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APPENDIX B 



In order to calculate the low-energy phase shifts we follow the procedure of approximating 
the interior wave function by a small fc^ expansion about its zero-energy solution. Thus, 
(11071) has the form 



rf,2+/±W)P±i = -fcV±i. (Bl) 



At a zero energy 



From fIBll) and flB2l) get, after integrating. 



0. (B2) 



Po,±i (^)pli {r) - Pi 1 (r)pQ ^1 (r) = -A;2 ! dsp^ipo^i , (B3) 

2 2 2 > 2 J 2 



since p,! ~ r^^+l provided r^a = 0, which we have assumed. Hence 

2 r^O r^O 



P±i , 



7^1 - (fci?)2r±i + 0{kR)\ (B6) 



Here 



V9,.po 



. 1 



7±i 



5 



where 7^1, r_|_i are denoted by 7^, F^, respectively, in Sees. VI and VII with a denoting 
j, M, ±i. When M = —j — i, Pq _i is given by fll53p and when M = j + ^, 



(B4) 



P±l ) R \ Po,±i P±i(^)Po,±i(^) 
Since /c^ is an analytic perturbation of p^i in (IB II) . make the expansion 

P±i W = Po,±i + P2,±i(0^' + 0{k'), (B5) 

and substitute this in (lB4p to get the interior logarithmic derivative 

rc}r.p_|_] 



Po„ 
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Calculation of the low-energy phase shifts then proceeds by equating the interior deriva- 
tive f lB6|) with the exterior derivative calculated from the small-A; expansion of 




obtained from (11131) and (11161) . 



42 



[1] G. V. Dunne, J. Hur, C. Lee and H. Min, Phys. Rev. Lett. D 94, 072001 (2005); Phys. Rev. 

D 71, 085019 (2005). 
[2] M. P. Fry, Phys. Rev. D 67, 065017 (2003). 

[3] M. P. Fry, Phys. Rev. D 62, 125007 (2000); D 72, 109903(E) (2005). 
[4] J. Schwinger, Phys. Rev. 128, 2425 (1962). 

[5] M. P. Fry, Int. J. Mod. Phys. A 17, 936 (2002). This reference contains numerous misprints. 

The original paper (hep-th/0111244) should be consulted. 
[6] E. Seller, in Gauge Theories: Fundamental Interactions and Rigorous Results, Proceedings 

of the International Summer School of Theoretical Physics, Poiana Brasov, Romania, 1981, 

edited by P. Dita, V. Georgescu, and R. Purice, Progress in Physics Vol. 5 (Birkhauser, Boston, 

1982), p.263. 

[7] E. Seller, Gauge Theories as a Problem of Gonstructive Field Theory and Statistical Physics, 

Lecture Notes in Physics Vol. 159 (Springer, Berlin, 1982). 
[8] E. Seller and B. Simon, Commun. Math. Phys. 45, 99 (1975). 

[9] B. Simon, Trace Ideals and their Applications, London Mathematical Society Lecture Notes 
Series 35 (Cambridge University Press, Cambridge, England, 1979). 
[10] B. Simon, Adv. Math. 24, 244 (1977). 
[11] S. L. Adler, Phys. Rev. D 16., 2943 (1977). 

[12] S. L. Adler, Phys. Rev. D 10, 2399 (1974); D 15, 1803(E) (1977). 

[13] I. C. Gohberg and M. G. Krem, Introduction to the Theory of Linear Nonselfadjoint Oper- 
ators, Translations of Mathematical Monographs Vol. 18 (American Mathematical Society, 
Providence, 1969), p. 163. 

[14] N. Dunford and J. T. Schwartz, Linear Operators, Part II (Interscience, New York, 1963). 

[15] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[16] R. Karplus and M. Neuman, Phys. Rev. 80, 380 (1950). 

[17] V. Weisskopf, Kong. Dans. Vid. Selsk., Math-fys. Medd. XIV No. 6 (1936). 

[18] M. P. Fry, Phys. Rev. D 55, 968 (1997). 

[19] S. L. Adler, Phys. Rev. D 6, 3445 (1972); D 7, 3821(E) (1973). 

[20] C. Itzykson, G. Parisi and J. B. Zuber, Phys. Rev. D 16, 996, (1977). 



43 



[21] R. Balian, C. Itzykson, J. B. Zuber and G. Parisi, Phys. Rev. D 17, 1041 (1978). 

[22] E. B. Bogomolny and Yu. A. Kubyshin, Sov. J. Nucl. Phys. 35, 114 (1982); Yad. Fiz. 35, 
202 (1982); reprinted in Large-Order Behaviour of Perturbation Theory, edited by J. C. Le 
Guillou and J. Zinn-Justin, Current Physics - Sources and Comments Vol. 7 (North Hohand, 
1990), p.425. 

[23] R. Musto, L. O'Raifeartaigh and A. Wipf, Phys. Lett. B 175, 433 (1986). 

[24] R. G. Newton, Scattering Theory of Waves and Particles (Springer, New York, 1982). 

[25] W. Heisenberg and H. Euler, Z. Phys. 98, 714 (1936). 



[26] G. V. Dunne, hep-th/0406216 



[27] L. S. Brown, R. D. Carhtz and C. Lee, Phys. Rev. D 16, 417 (1977). 
[28] G. 't Hooft, Phys. Rev. D 14, 3432 (1976); D 18, 2199(E) (1978). 



44 



